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ABSTRACT 

We present an X-ray absorption line spectroscopic study of the large-scale 
hot interstellar medium (HISM) in the Galaxy. We detect Ne IX K a absorp- 
tion lines in the Chandra grating spectra of seven Galactic low-mass X-ray bi- 
naries. Three of these sources also show absorption of VII K a , VII Kg, 
and/or VIII K a . Both centroid and width of the lines are consistent with a 
Galactic HISM origin of the absorption. By jointly fitting the multiple lines, 
accounting for line saturation and assuming the collisional ionization equilib- 
rium, we estimate the average absorbing gas temperature as ~ 2.4 ± 0.3 x 10 6 K 
(90% confidence errors). We further characterize the spatial density distribution 
of the gas as 6.4l^; 7 exp(— |z|/1.2ij;° kpc) x 10~ 3 cm -3 (a disk morphology) or 
[1 + (R/2.3ii;i kpc) 2 ]" 1 x 10" 3 cm -3 (a sphere morphology), where z and 
R are the distances from the Galactic plane and Galactic center (GC) respec- 
tively. Since nearly all the sight-lines with significant absorption lines detected 
are somewhat toward GC and at low Galactic latitudes, these results could be 
severely biased. More observations toward off-GC sight-lines and at high latitudes 
are urgently needed to further the study. Nevertheless, the results demonstrate 
the excellent potential of X-ray absorption line spectroscopy in the study of the 
HISM. 

Subject headings: Galaxy: halo — X-rays: ISM — X-rays: individual (NGC 6624, 
GX 9+9, Ser X-l, Cyg X-2, V801 Ara, V926 Sco, GX 349+2) 



1. Introduction 

It has long been theorized that a major, possibly dominant, phase of the interstellar 
medium (ISM) in both the disk and the halo of the Galaxy is gaseous at temperature T ~ 10 6 
K (Spitzer 1956; Cox & Smith 1974; McKee & Ostriker 1977; Heiles 1987; Ferriere 1998). 
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The presence of this rarefied hot ISM (HISM) component affects the geometry and dynamics 
of the cool phases of the ISM, the propagation of cosmic rays and UV/soft X-ray photons, 
the Galactic disk-halo interaction, the distribution of metal abundances, and the feedback 
to the intergalactic medium. There is, however, no consensus on how much hot gas exists, 
how it is distributed in the Galaxy, and what thermal, chemical, and ionization states it is 
in. 

Hot gas with temperatures of T > 10 6 K may be traced by X-rays. Indeed, it has 
been studied extensively with broad-band X-ray observations, producing all-sky maps of the 
diffuse soft X-ray background (SXB, e.g., Snowden et al. 1997). Furthermore, emission lines 
such as C VI, O VII, and O VIII have been detected in a high spectral resolution observation 
made with microcalorimeters flown on a sounding rocket, confirming the thermal origin for 
much of the background (McCammon et al. 2002). However, the X-ray emission carries 
little distance information, and its interpretation is typically subject to large uncertainties 
in line-of-sight absorption. Therefore, it is very difficult, if not impossible, to infer the spatial 
and physical properties of the hot gas from X-ray emission measurements alone. 

Alternatively, one can study the HISM by observing its absorption against bright back- 
ground X-ray sources. This capability is now provided by the grating instruments aboard 
Chandra and XMM-Newton X-ray Observatories (Weisskopf et al. 2000; Jansen et al. 2001). 
Indeed, the highly-ionized oxygen/neon absorption lines (O VII, O VIII, and/or Ne IX), 
consistent with no velocity shift, have been detected in the spectra of several bright active 
galactic nuclei (AGNs): PKS 2155-304 (Nicastro et al. 2002; Fang et al. 2002), 3C 273 
(Fang et al. 2003), and MKN 421 (Nicastro 2003; Rasmussen et al. 2003). Wang et 
al. (2005) have further detected narrow O VII and Ne IX absorption lines in the spectrum 
of LMC X-3; the equivalent widths (EWs) of these lines are similar to those seen in the AGN 
spectra, suggesting that the bulk of the absorbing material is within the 50 kpc distance of 
LMC X-3. In addition, Futamoto et al. (2004) have detected O VII and O VIII absorption 
lines in the sight-line toward a Galactic low mass X-ray binary (LMXB) X1820-303; the 
equivalent widths of the lines are substantially higher than those seen in the AGN spectra, 
consistent with a stronger diffuse X-ray emission observed in the Galactic bulge region, where 
this LMXB is located. 

The detection of highly- ionized species in absorption lines potentially provides a powerful 
tool in the study of the HISM. For ease of reference, Fig. 1 presents the relative ionization 
fraction of oxygen and neon as a function of temperature, for a plasma in a collisional 
ionization equilibrium (CIE) state. The X-ray absorption lines are sensitive to the gas over 
the entire temperature range expected for the HISM. 

• An individual absorption line alone, if resolved, gives a direct measurement of the 
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Fig. 1. — The ionization fractions of oxygen (thin lines) and neon (relative to all oxygen, 
thick lines) as a function of the CIE plasma temperature (Arnaud & Rothenflug 1985). 

kinematics and column density of the X-ray-absorbing gas, independent of the filling 
factor of the gas along the line of sight. 

• Two or more absorption lines from the different transitions of the same ion, i.e., K a , 
Kg, etc., even if not resolved, may also be used to constrain the kinematics and column 
density of the gas. 

• Two or more lines from the same element (s) but from different ionization states may 
provide diagnostics of the thermal and ionization states of the gas. 

• Multiple absorption lines of different elements (and transitions) may further allow to 
measure their relative abundances. 

• If absorption lines are detected along many sight-lines and if the distances to the 
background sources are known, one may then characterize the spatial distribution of 
the gas. 

• Because the absorption is relative to the local continuum level, all such absorption line 
measurements are independent of line-of-sight cool gas absorption! 

• Of course, a joint analysis of the X-ray absorption line(s) with (line and/or broad-band) 
emission measurements may further give additional constraints on such important pa- 
rameters as the average volume density and the filling factor of hot gas. 

We have carried out a systematic Chandra archival study of the X-ray absorption lines in 
the spectra of Galactic LMXBs. Because of their brightness and intrinsic spectral simplicity 
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(relatively free from the complexity which could be caused by the stellar wind of a massive 
companion; e.g., Schulz & Brandt 2002), LMXBs are excellent background sources for the 
study of the intervening X-ray-absorbing gas. We find that existing archival observations 
are already quite useful, although they were typically not optimized for detecting absorption 
lines, in terms of both the instrument setup and the exposure time. In addition to the work 
by Futamoto et al. (2004), there have been a few reported X-ray absorption line studies 
that are based on the grating observations of LMXBs. However, these studies are focused on 
the absorption lines produced either by the cool phase of the ISM (e.g., Paerels et al. 2001; 
Juett et al. 2004) or by gas intrinsic to the binaries (e.g., Miller et al. 2004). 

Here we report results from our archival study. We first describe our selection of the 
LMXBs and the Chandra observations in §2 and then introduce our multiplicative absorption 
line model in §3. In §4, we present the results of our absorption line measurements. We 
discuss the origin (§5) and spatial distribution (§6) of the X-ray-absorbing gas and make 
comparisons with other relevant observations (e.g., O VI and diffuse X-ray emission; §7). 
Finally in §8. we summarize our results and conclusions. Throughout the paper, errors are 
quoted at the 90% confidence level. 

2. Source Selection and Observations 

Our study use 17 Chandra grating observations of 10 LMXBs, available in the archive 
in April 2004 (Table 1; Fig. 2). These sources are selected with two criteria: 1) a Galactic 
latitude |6| > 2° to minimize the effect of soft X-ray absorption by cool gas, and 2) a high 
signal-to-noise ratio > 7 per bin at ~ 0.9 keV to have a good probability for a positive 
absorption line detection. We exclude those sources with absorption or emission lines that 
have been identified as intrinsic to the binaries (e.g., GX339-4, Miller et al. 2004; Her X-l, 
Jimenez-Garate et al. 2002). 

Chandra carries two high spectral resolution instruments: the high energy transmission 
grating (HETG, Markert et al. 1995; Canizares et al. 2000) and the low energy transmission 
grating (LETG, Pease et al. 2002). The HETG consists of two grating assemblies, the 
high energy grating (HEG) and the medium energy grating (MEG). The LETG can be 
operated with either the advanced CCD imaging spectrometer (ACIS) or the high resolution 
camera (HRC), whereas the HETG works only with the ACIS. The energy resolution of 
the MEG, HEG, and LETG is 0.023, 0.012, and 0.05 A (FWHM), respectively 1 . For 



1 For detailed Chandra instrumental information, please refer to Proposers' Observatory Guide 
http://cxc.harvard.edu/proposer/POG/index.html 
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Table 1. Source Parameters 







(I, b) 


D 


z 


P 


A 




Source 


Other Name 


(cleg) 


(kpc) 


(kpc) 


(hr) 


(R ) 


Ref. b 


X1820-303 


NGC 6624 


(2.79, -7.91) 


7.6 


1.0 


0.19 


< 0.1 


1,2,3 


X1728-169 


GX 9+9 


(8.51, 9.04) 


5.0 


0.8 


4.2 


1.6 


1,4,5 


X1837+049 


Ser X-l 


(36.12, 4.84) 


8.4 


0.7 


13 


~3.4 a 


1,4 


X2142+380 


Cyg X-2 


(87.33, -11.22) 


7.2 


1.4 


235 


26.9 


1,9,15 


XI 118+480 




(157.66, 62.32) 


1.8 


1.6 


4.1 




10, 11 


X0614+091 




(200.88, -3.36) 


< 3 


< 0.2 


< 1 




12 


X0929-314 




(260.11, 14.22) 


> 5 


> 1-2 


0.72 




13,16 


X1636-536 


V801 Ara 


(332.92, -4.82) 


6.5 


0.5 


3.8 


1.6 + 0.3 


1,4,7 


XI 735-444 


V926 Sco 


(346.05, -6.99) 


7.1 


0.9 


4.6 


1.75 + 0.02 


1,4,8 


XI 702-363 


GX 349+2 


(349.10, 2.75) 


5.0 


0.2 


22 


~4.8 a 


14, 4 



Note. - D, z, P, and A are the distance, vertical distance away from the Galactic 
plane, orbit period, binary separation of the LMXBs. 

a Estimated in this work by assuming that the companion is a main sequence star with a 
mass of M c ~ 0.4 M , overflowing its Roche lobe and being accreted by a compact object 
with a mass of M x ~ 1.4 M (Paczynski 1971). 

References: Van Paradijs 1993; 2 Kuulkers et al. 2003; 3 Stella et al. 1987; Christian 
& Swank 1997; 5 Schaefer 1990; 6 Homer et al. 1996; 7 Lawrence et al. 1983; 8 Smale & 
Corbet 1991; 9 Cowley et al. 1979; 10 Cook et al. 2000; n McClintock et al. 2001a; 12 Brandt 
et al. 1992; 13 Markwardt et al. 2002; 14 Wachter k Margon 1996; 15 Vrtilek et al. 2003; 
16 Galloway et al. 2002; 
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Fig. 2. — Directions of X-ray sources in the present study, shown on a ROSAT all-sky 
diffuse 3/4-keV-band intensity map (in the Aitoff projection; Snowden et al. 1997). These 
sources are the 10 LMXBs (filled circle) listed in Table 1, plus the two extragalactic sources, 
MRK 421 and LMC X-3 (filled diamond). The blue symbols represent those with significant 
absorption line detections, whereas the white symbols mark ones without a detection. 

ease of reference, we show in Fig. 3 the effective area of the telescope/instruments as a 
function of photon wavelength. With the expected range of hot gas column density, the 
most promising absorption lines are at the wavelengths of ~ 33 A (C VI), 20 A (O VII and 
O VIII), and 13 A (Ne IX and Ne X). The rest-frame line wavelengths are marked in Fig. 3. 
Note that the effective area of the LETG/ACIS is small at ~ 33 A, where the cool gas 
absorption of the continuum is also severe. Furthermore, for a LETG/HRC spectrum, the 
higher order confusion could be severe at this wavelength and sorting orders may also cause 
some potential uncertainties. Therefore, the C VI is more difficult to detect than the other 
lines. The LETG/ACIS and LETG/HRC, with an effective area of ~ 20/30 cm 2 at 20/13 A, 
should typically be the optimal combination for detecting O VII, O VIII, and Ne IX lines 
in one observation. The actual sensitivity of the LETG/ACIS also depends on the offset 
pointing, which determines where the source is placed on the detector and whether the 
lines are detected in the front-illuminated CCDs or the back-illuminated CCDs. For HETG 
observations, only the data from the MEG are considered here because it is substantially 
more sensitive than the HEG. Also because of the higher resolution of the MEG (by a factor 
of ~ 2), it can be more sensitive for detecting narrow Ne IX and Ne X lines than the LETG. 
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Table 2. Chandra Observations 







Obs. Date 




Exposure 


Source 


ObsID 


/ /ll/ \ 

(mm/dd/yyj 


Grating/Detector 


(ks) 


XI 820-303 


98 


03/10/00 


LETG/HRC 

J_J 1 J _L V.I / -L J. J. I V. 


15.12 




1021 


07/21/01 


HETG/ACIS 


9.70 




1022 


09/12/01 


HETG/ACIS 


10.89 


X1728-169 


703 


08/22/00 


HETG/ACIS 


22.44 


X1837+049 


700 


06/13/00 


HETG/ACIS 


78.06 


X2142+380 


87 


04/24/00 


LETG/HRC 


30.16 




111 


11/11/99 


LETG/ACIS 


6.2 




1016 


08/12/01 


HETG/ACIS 


15.13 




1102 


09/23/99 


HETG/ACIS 


28.98 


X1118+480 


1701 


04/18/00 


LETG/ACIS 


27.83 


X0614+091 


100 


11/28/99 


LETG/HRC 


26.23 


X0929-314 


3661 


05/15/02 


LETG/ACIS 


17.96 


X1636-536 


105 


10/20/99 


HETG/ACIS 


29.78 




1939 


03/28/01 


HETG/ACIS 


27.06 


XI 735-444 


704 


06/09/00 


HETG/ACIS 


24.91 


X1702-363 


715 


03/27/00 


HETG/ACIS 


10.98 




3354 


04/09/02 


HETG/ACIS 


35.21 
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Among the 17 observations, only six (on five sources) used the LETG (Table 2). All of 
them, except for the short exposure on X2142+380 (ObsID 111, a calibration observation), 
have been analyzed for detecting absorption lines in previous studies. No significant highly 
ionized oxygen or neon absorption line is detected for X1118+480 (McClintock et al. 2001b), 
X0614+091 (Paerels et al. 2001), and X0929-314 (Juett et al. 2003). Our re-analysis of 
the data confirms this result. The detection of O VII K a and Kp, O VIII K a , and Ne IX K a 
absorption lines has been reported for X1820-303, and an upper limit to the EW of the O VII 
absorption line has been set from the LETG spectrum of X2 142+380 (Futamoto et al. 2004). 
Our analysis of these two sources used both LETG and HETG observations to achieve higher 
sensitivities. We report here for the first time the detection of the highly ionized absorption 
lines toward X1728-169, X1837+049, X1636-546, X1735-444, and X1702-363. 

We re-process the observations, using the standard CIAO software (version 3.02) with 
the calibration database CALDB (version 2.25) to extract source and background spectra. 
For each observation we calculate the auxiliary response functions (ARFs) by running the 
CIAO thread fullgarf for the positive and the negative grating orders, and adopt the response 
matrix files (RMFs) from the CALDB directly. All the HETG/ACIS observations had the 
zeroth order images either severely piled-up or intentionally blocked to avoid telemetry sat- 
uration. We use the intersection of the two grating arms (both the HEG and the MEG) and 
the readout "streak" to determine the position of the source in the detector and to calibrate 
the wavelength of the grating spectra (e.g., Schulz et al. 2002). We then co-add the positive 
and the negative order spectra to improve the counting statistics. We further combine the 
spectra from multiple observations of a source, using the CIAO thread add_grating_spectra. 



3. Implementation of an Absorption Line Model 

Published X-ray absorption line analyses typically use the following procedure (e.g., 
Futamoto et al. 2004): (1) adopting a smooth function (e.g., a power- law) to characterize 
the spectral continuum and an additive negative Gaussian function to account for the ab- 
sorption line profile; (2) fitting the observed spectrum to constrain the parameters of the 
functions; (3) using the fitted EW of the line in a curve-of-growth analysis to estimate the 
ionic column density. This procedure is simple but typically only adequate for individual 
unsaturated line analysis. The procedure also does not use all the information available in 
the observed spectrum: details of the observed line shape and/or the intrinsic connections 
between multiple lines. Although a joint analysis of the ionic column densities inferred from 
the analysis of individual lines can, in principle, be used to constrain the related physical 
parameters such as temperature and ionic abundances, it is generally difficult to correctly 
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propagate the non-Gaussian errors of the estimated parameters. 

We have implemented a general model (referred here as absline) for absorption line 
fitting in the X-ray spectral analysis software package XSPEC. Following the description of 
the absorption process given by Rybicki & Lightman (1979) and Nicastro et al. (1999), the 
radiation transfer at photon energy e can be expressed as 

1(e) = I c {e)e-^\ (1) 

and 

7TP 2 

r( £ )=iV J f lu <j>(e), (2) 

m e c 

where m e is the electron mass, c is the speed of light, and is the oscillator strength of 
the electron transition from a lower level to an upper level. The ionic column density iVj is 
a function of the reference element (hydrogen as the default) column density N H and the 
plasma temperature T, 

N* = N H fJ c (T), (3) 

where f a is the corresponding element abundance relative to the reference element and fc(T) 
is the ionic fraction as a function of T (assuming that the absorbing gas is in a CIE state). 
The normalized Voigt function <j)(e) is a convolution of the intrinsic Lorentz line profile with 
a Doppler broadening (assumed to be Gaussian): 

^ = L tfF^Tw ^ D d(Ac > (4 > 

' H(a,u), (5) 



if (a, u) 

where 



a f + °° e~y 

- / t — ^2~i — 2 d y ' 6 

J-oo [u - vY + a 2 



u = ( 7 ) 



4nAe D Ae D 

Here 7 is the natural broadening damping factor, e\ is the systemic (rest frame if redshift 
z=0) energy of the line, and Aeo is the Doppler width 

Ae D = (8) 

in which 

b v = (—+e) 1/2 , o) 
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where rrii is the ionic mass, £ is the velocity dispersion due to extra-broadening (e.g., tur- 
bulence). When u is small (in the core region of a line), H(a,u) ~ e - " 2 gives a Gaussian 
profile, whereas when u is large (in the wings of a line) if (a, u) ~ a/nu 2 is close to a Lorentz 
profile. 

We define e~ T ( £ \ part of Eq. 1, as our multiplicative absline model, which is specified 
by five parameters: the plasma temperature T, central line energy ion velocity dispersion 
b v , the hot-phase hydrogen column density Nh (or JVj for a chosen ion), and the metal 
abundance f a . Other parameters (fi u and 7) are given for a specific line (Table 3). While 
ei is very sensitive to the line centroid in a fit, b v and iVj can be constrained by the shape 
and intensity of the line. When multiple lines are present, one can conduct a joint fit, which 
may also allow for the estimation of T, Nh, and/or f a and their uncertainties. 

Note that in the absline model, b v is used as an independent parameter (not constrained 
by Eq. 9), and can be smaller than the CIE thermal broadening (2/cT/mj) 1 / 2 . Therefore, 
with the adopted Voigt function, the absline model can be applied as well to a photo-ionized 
plasma where b v ~ 0. The inclusion of the Lorentz profile in calculating the line profile is 
important when the Doppler broadening is small (e.g., in an over-cooled or photo-ionized 
plasma) and/or when the column density is large. Using the Lorentz profile alone gives a 
firm upper limit to the ionic column density (e.g., Futamoto et al. 2004). 

Fig. 4 illustrates the difference between the additive Gaussian and the absline models. 
The additive Gaussian model is a good approximation to the absline model only when r is 
small (< 1). When To is large (> 1), the additive Gaussian profile deviates significantly from 
the absline profile and the inferred velocity dispersion be, for example, can be substantially 
overestimated. 



4. Analysis and Results 

Our absorption line analysis is based on the data in 12-14, 18-20, and 20-22 A band- 
passes, which embrace the Ne X K a and Ne IX K a , O VII and O VIII K a , and O VII K Q 
lines, respectively (Table 3). The chosen band width (2 A) is a compromise between maximiz- 
ing the counting statistics and reducing the effect of the potential deviation of the spectral 
continuum from our power-law characterization. The subtracted background contributes 
< 5% to the source counts in the wavelength bandpasses for all sources. We apply both 
the commonly used additive Gaussian model and the new absline model to the absorption 
line analysis of the 10 selected sources (Tables 1-2). The models, including both the con- 
tinuum and the line(s), are convolved with the instrument responses in XSPEC (version 
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Fig. 3. — The effective area of the first-order spectra of the HEG (dashed line), the MEG 
(solid line), the LETG/ACIS (dash-dotted line), and the LETG/HRC (dotted line; with 
unsorted orders). The vertical lines mark the rest frame wavelengths of highly ionized neon, 
oxygen, and carbon species considered in this work. The ACIS data were adopted from 
Chandra cycle 4. 



Table 3. Absorption Line Parameters 





rest \/si 




7 


Line 


(A/eV) 


flu 


(10- 2 eV) 


Ne X K Q 


12.134/1021.79 


0.416 


2.61 


Ne IX K a 


13.448/921.95 


0.657 


3.35 


VII K p 


18.629/665.55 


0.146 


0.39 


VIII K a 


18.967/653.69 


0.277 


1.06 


VII K a 


21.602/573.95 


0.696 


1.37 



Note. - The f u is the transition oscillator 
strength and 7 is the natural broadening damping 
factor. The values of Ne IX are adopted from Be- 
har & Netzer (2002) and the others from Verner et 
al. (1996). 
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Fig. 4. — Comparison between the absline model of a Ne IX K a absorption line (solid line) 
and the best-fit additive Gaussian model (dotted). The absline model assumes a velocity 
dispersion 6^=45 km s -1 (thermal velocity at T~ 2.4 x 10 6 K) as well as log(N NcIX )=15.0, 
15.5, 16.0, and 16.5 cm -2 (for the panels from left to right), corresponding to optical depths 
of r =0.29, 0.91, 2.89, and 9.10 at the line centroid (921.95 eV ~ 13.448 A). The b v = 
Lorentz profiles (dashed) are also shown for reference in the panels. The limitations of the 
additive Gaussian model are characterized by the best fit b G values as they significantly 
differ from 6^=45 km s -1 . 

11.3.1) before they were fitted to the observations. We choose the signal-to-continuum-noise 
ratio S/N = 3a (corresponding to a false detection probability ~ 0.1%) as the detection 
threshold. Seven sources show significant Ne IX K a absorption lines in the MEG spectra; 
three of these sources also show O VII Kg, O VIII K Q , and/or O VII K Q absorption line(s). 
The O VII K a absorption line is not significant in the LETG/HRC spectrum of X2142+380, 
which is consistent with the previous work by Futamoto et al. (2004), but we detect sig- 
nificant Ne IX K a and O VII Kg absorption lines in the more sensitive MEG spectra. All 
the detected lines are unresolved at 90% confidence level. The lack of significant absorption 
lines in the spectra of X1118+480, X0614+091, and X0929-314 is consistent with the results 
from the previous studies as mentioned in §2. No Ne X K a absorption line is detected in 
any of the 10 sources. 

For the detection of individual lines, it is convenient as a first pass to use the Gaussian 
model. We estimate the line centroid energy and EW for each detected line. Fig. 5 shows the 
Gaussian model fits to the Ne IX K a lines. For each sight-line, we also estimate an upper 
limit to the EW of the Ne X K a line by fixing its line centroid energy at the rest frame 
energy and jointly fitting its width with the Ne IX K a line. 
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Fig. 5. — The Ne IX K a absorption lines in the spectra of the 10 sources. The solid line 
in each panel represents a Gaussian plus power-law model folded with the instrumental 
response; the vertical dashed line marks the rest frame wavelength of the line. The apparent 
broad line features in (d), (e), and (g), are partly due to dip in the instrument effective 
area near the line central wavelength; the short vertical arrows indicate other instrumental 
features. The bin size is 12.5 mA for LETG observations (panels e, f, g) and is 5.0 mA for 
HETG observations (the other panels). 
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Using the absline model, we jointly fit multiple lines of each source. Table 4 marks the 
included lines. In such a fit, we link the model parameters b v , T, and Nh of the lines; the 
slight dependence of b v on the different ion mass is neglected (< 5 km s -1 ; see Eq. 9). Because 
of the limited number and quality of the line detections, we use a fixed ISM abundance 
(Wilms et al. 2000) and assume the CIE for the X-ray-absorbing gas. For those lines that 
are not detected, we fix ei to their rest-frame energies. Fig. 6 presents the velocity profiles 
of the jointly-fitted multiple lines. 

Our results are summarized in Table 4. The sources X1118+480, X0614+091, and 
X0929-314 are not included in the table; their poor data quality (in terms of both counting 
statistics and spectral resolution) does not give any meaningful upper limits to either EWs 
or the ionic column densities. The S\ values and their uncertainties constrained from absline 
model are consistent with those obtained from the Gaussian model and are thus not repeat- 
edly presented in the table. In addition to setting the upper limits to b v from the observed 
lines, the joint-fits further constrain the lower bounds on b v for four of the sources (Table 4). 
If b v were smaller, the (unresolved) line would then be deeper in order to provide the same 
observed absorption. The degree of saturation To varies among the jointly-fitted lines and 
therefore affects their relative strengths, which sets a lower bound on b v . Table 4 further 
includes a hot gas hydrogen density %, estimated from assuming a uniform distribution 
with a unity filling factor along each sight-line. 

Considering the potential uncertainties in the systematics of the spectral resolution 
calibration, we have also estimated the Ne IX column density N^ cIX by fitting individual 
Ne IX K Q lines with the natural broadening only (§3). This "firm" upper limit to the Ne IX 
column density is presented in the last column of Table 4. 

The joint-fits give direct estimates of Nh and T, although the uncertainties are large 
along some sight-lines. We define a mean measurement p and its 90% upper and lower 
uncertainties Ap ± of a parameter as 



where pi, Ap +i and Ap_j are the measurement and its corresponding 90% upper and lower 
uncertainties along the ith sight-line with the absorption line(s) detected. We obtain T = 
2.4(2.1, 2.7) x 10 6 K, N H = 8(6, 10) x 10 19 cm" 2 , and n H = 3.6(3.0, 4.5) x 10~ 3 cm" 3 . 

The plasma properties along the sight-line of X1820-303 are mostly consistent with those 
obtained by Futamoto et al. (2004) from a curve-of-growth analysis. However, they prefer 
a large b v (> 200 km s -1 ), arguing that otherwise the inferred N H would be comparable 
to, or even larger than, the neutral hydrogen column density observed in the field. Our 
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Fig. 6. — The velocity profiles of the joint-fitted absorption lines. The transition and centroid 
wavelength of each line as well as the used instrument are listed in each panel. The solid 
histogram of each panel represents the best-fit absline-modi&ed power-law model folded with 
the instrumental response. The vertical dashed line marks the zero velocity. A negative 
velocity indicates a blue shift. The bin sizes are the same as in Fig. 5. 
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Fig. 6. — Continued. 



Table 4: Results of Gaussian and absline model fits 



S/N 

Source Line (<j) 


Gaussian model 
cz b G EW 
(km/s) (km/s) (cV) 


absline model 

K log(T) log(N H ) n H 1 °g( N & e rx) 
Joint (km/s) (K) (cm -2 ) (10~ 3 cm~ 3 ) (cm" 5 ) 


X1820-303 Ne X K a < 1 
Ne IX K a 8.6 
O VII Kp < 1 
O VIII K a 6.0 
O VII K Q 10.6 


< 0.09 

-89(-112, 89) <346 0.43(0.28, 0.63) 

-86(-201, 152) <490 0.65(0.34, 0.98) 
175(-120, 125) <347 0.58(0.41, 0.79) 


V 
V 

191(62, 346) 6.4(6.2, 6.5) 20.0(19.8, 20.2) 4.3(2.7, 6.8) <18.54 


X1728-169 Ne X K Q < 1 
Ne IX K a 4.9 
O VII Kp 2.0 
O VIII K a 2.1 


< 0.36 

-134(-223, 134) <532 0.30(0.14, 0.52) 
-261(-342, 264) <375 0.39(0.14, 0.65) 
178(-298, 617) <531 0.43(0.16, 0.70) 


V 

< 375 6.3(6.0, 6.5) 19.9(19.7, 22.2) 5.1(3.2, 1027) <18.37 


X1837+049 Ne X K Q < 1 
Ne IX K a 9.9 
O VII Kp < 1 
O VIII K a < 1 


< 0.11 

22(-89, 89) <366 0.42(0.37, 0.47) 


V 
V 
V 

■J 

280(111.443) 6.45(6.37, 6.54) 20.1(19.9, 20.2) 4.8(3.1, 6.1) <18.46 


X2142+380 Ne X K Q < 1 
Ne IX K Q 4.8 
O VII Kp 4.9 
O VIII K a < 1 
O VII K a < 1 


< 0.11 

-223(-267, 290) <797 0.25(0.14, 0.40) 
l(-263, 285) <369 0.24(0.01, 0.41) 


V 
V 

•J 
-J 

383(30, 664) 6.3(6.2, 6.4) 19.7(19.5, 19.9) 2.3(1.4. 3.6) <18^) 


X1636-536 Ne X K a < 1 
Ne IX K a 8.4 
O VII Kp < 2 
O VIII K a < 2 


< 0.09 

-67(-44, 45) <213 0.38(0.27, 0.49) 


i 
•j 

<105 6.3(5.7, 6.8) 20.2(19.5, 22.3) 7.9(1.6, 995) <18.28 


X1735-444 Ne X K a < 1 
Ne IX K a 9.0 
O VII Kp < 2 
O VIII K a 3.8 


< 0.31 

-22(-156, 200) <724 0.78(0.44, 0.84) 
-168(-611, 375) <746 <0.90 


V 

•J 
•J 
•J 

369(58, 638) 6.5(6.3, 6.6) 20.2(20.0, 20.4) 7.2(4.6. 11) <18.84 


X1702-363 Ne X K a < 1 
Ne IX K a 5.5 
O VIII Kp < 1 
O VIII K a < 1 


< 0.28 

-201(-289, 335) <593 0.50(0.17, 0.60) 


V 

•J 

<899 6.0(5.6, 7.1) 19.9(19.4, 20.3) 5.1(1.6, 13) <18.62 



Note. — The VII K Q line is only searched in the spectra of LETG observations (Table 2). 
A negative cz value indicates a blue shift. In a Gaussian model fit, be is \/2 times the standard 
deviation (mimicking the velocity dispersion b v in an absline model fit). The "Joint" column marks 
those absorption lines utilized in the joint-fits with the absline model. The nn is the averaged hot 
gas density derived from the column density Nh divided by the source distance D listed in Table 1. 
All limits are at the 90% confidence level. See the text for details. 
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b v value can be as low as 62 km s _1 , and the hot gas column density is well below the 
neutral hydrogen column density of 1.9 x 10 21 cm -2 (Bohlin, Savage, & Drake 1978). This 
discrepancy is probably due to a significant deviation of their used additive Gaussian model 
from the proper multiplicative absline profile (Fig. 4) and to the inaccurate error propagation 
of the temperature in their analysis (§3). 

5. Origin of the Hot Gas 

In principle, the X-ray-absorbing hot gas could be local to the binary systems. To 
account for the large column densities estimated above, a plausible scenario might be the 
photo-ionized winds from the accretion disks that are presumably responsible for the X-ray 
continuous emission. Such a scenario, if confirmed, would be interesting in its own right. 
However, we find that it has serious difficulties: 

• An absorption line produced in an accretion disk wind should be blue-shifted and 
broadened with a magnitude comparable to the escape speed of the accreting compact 
object (e.g., Ueda et al. 2004). The line shift, for example, can be characterized by 
AA/A = v w /c ~ (2r 9 /r i ) 1 / 2 , where v w is the wind velocity, r ; is the starting radius 
where the absorption line is produced, r g = GM lc 2 is the gravitational radius. Since 
all the detected X-ray absorption lines are consistent with cz ~ 0, using the wavelength 
accuracy 0.01 1A of MEG as AA and assuming the compact object mass M = 1.4M Q , 
we obtain ri > 9 R Q . This required r ; value is much greater than the binary separation 
A (Table 1) for all the sources with absorption lines detected except for X2 142+380. 

• The Ne IX absorption should arise only in a region with a proper ionization state. If 
the electron density in the wind can be approximated as n(r) oc l/r a (a > 1), the 
ionization parameter can then be written as, 

U = L x / [n(r)r 2 ] = L x r a ~ 2 / [N H (a - l)^" 1 ] , 

where r > r\. According to Kallman & McCray (1982), Ne IX is the dominant 
ionization state of neon only when log(C/)< 2 in both optical thin and thick cases. We 
estimate the luminosity of the individual sources to be in the range of 2 — 7x 10 37 ergs s _1 
(assuming the source distances in Table 1). Further assuming a Ne IX fraction ~ 0.5 
and taking N H ~ 10 23 cm -2 (corresponding to the maximum value of the N^ eIX ; 
Table 4), a — 2 (constant wind velocity), we obtain ri > 7.5 R & . Again, the required 
ri is much larger than the binary separation A of all the sources except for X2 142+380. 
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• We should also expect additional signatures of the wind. If the wind is launched at 
a radius r w away from the compact object, the hot hydrogen column density between 
r w and ri is then N w = Nh [{i r i/r w ) a ~ 1 — 1], where Nh is inferred from the X-ray 
absorption lines. Taking r w = 2 x 10 8 cm (corresponding to ~ 10 3 r g ,3 x 10~ 3 i? Q ; Li 
& Wang 1999), r\ ~ lOi? , and a = 2, we got N w ~ 10 3 iV#. Such a wind should then 
result in strong emission lines, which are, however, absent in the spectra of the sources. 
Furthermore, because the ionization state should increase with decreasing radius, we 
should observe a column density of Ne X comparable to, or greater than, that of Ne IX. 
But we do not detect any significant Ne X absorption line in any of the spectra. 

Based on these arguments, even though not conclusive (there are other possibilities such 
as a clumpy wind), we find that the X-ray line absorption is unlikely to be associated with 
the LMXBs. An interstellar origin of the absorption is most likely and is assumed in the 
following discussion. 



6. Spatial Distribution of the Hot Gas 

To tighten the constraints on the distribution of the X-ray-absorbing gas, we include the 
results from the absorption line study of the X-ray binary LMC X-3 (Wang et al. 2005) and 
an AGN MRK 421 (Yao et al. 2005). The Nh measurements, as listed in Table 5, are made 
in the same way as the 7 Galactic LMXBs reported above. With this set of 9 column density 
measurements, though still quite sparse, we attempt to characterize the overall spatial scale 
of the X-ray-absorbing gas, by assuming two extreme distribution geometry, an infinite disk 
or a sphere. 

For the disk distribution, we assume 

n(z) =n e-W /z \ (10) 

where no is the mean gas density at the Galactic plane and Zh is the vertical scale height. 
The column density can then be expressed as 

A^ = n ^(l-e-l 2 l/^/sin|&|. (11) 

A x 2 fit to the data (Fig. 7) gives n = 6.4(4.7,8.8) x lO^cm" 3 and ^=1.2(0.7, 2.2) kpc 
with the best-fit x 2 /d.o.f. = 7.2/7. Here, we have set the \z\ value of MRK 421 to be ~100 
kpc for ease of fitting and plotting; the exact value is irrelevant as long as it is much larger 
than Zh- Furthermore, the 90% confidence uncertainties of the individual Nh measurements 
have been scaled by a factor of 0.61 to mimic the required 1 a errors in the x 2 fit- 



Table 5. Results of the Sight-lines toward Extragalactic Sources 





(l,b) 


D 


log(N H ) 


Source 


(cleg) 


(Mpc) 


(cm -2 ) 


LMC X-3 


(273.58, -32.08) 


0.05 


19.7(19.2, 20.5) 


MRK 421 


(179.83, 65.03) 


118 


19.4(19.3, 19.6) 



Note. - The Nh values were adopted from Wang 
et al. (2005) and Yao et al. (2005) respectively. 
The distance to MRK 421 was calculated by assuming 
H = 75 km s -1 Mpc -1 . 
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Fig. 7. — Left panel: The vertical hot gas column density as a function of the distance from 
the Galactic plane. Diamonds mark the measurements toward the LMXBs and triangles 
represent the line-of-sights to the extragalactic sources (see the text for details). The dashed 
line is the best-fit exponential function with n = 6.4 x 10~ 3 cm~ 3 and z h = 1.2 kpc (Eq. 11). 
Right panel: Zh vs. n confidence contours at the 68%, 90%, and 99% confidence levels. 
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For the spherical distribution, we assume 

n(R)=n c [l + (R/R h ) 2 }-\ (12) 

adopted from the so-called /5-model n(R) = n c [1 + (R/ Rh) 2 ]~ 3f3 ^ 2 with /3 = 2/3 (Sarazin 
1988; Jones & Forman 1984), where n c is the hot gas density at the Galactic center (GC) 
and Rh is the scale radius. The hot gas column density toward a source with Galactic 
coordinates (/, b) and a distance D (Fig. 8) can be calculated via the following integration, 

-l 

dr (13) 




where a 2 = 1 + Rq sin 2 9/ R 2 h) x = (r — R cos6)/R h , R is the distance between the Sun and 
GC (taken as 8 kpc in this work), and 9 is the angular separation between GC and the source 
(cos# = cos / cos b). The fit (x 2 /d.o.f = 9 .4/7) of the data with this model (Fig. 9) is slightly 
worse than that with the disk model. The best-fit parameters are n c = 6.2(4.1, 10.0) x 10~ 3 
cm' 3 , R h = 2.3(1.2,3.9) kpc. 

With the fitted n and z h , or n c and R h , we estimate the total hot gas mass as ~ 2.7 x 10 s 
M Q and ~ 1.2 x 10 8 M for the disk and the sphere characterizations within a 15 kpc radius. 
This mass estimate is very uncertain as it depends sensitively on the assumed characteriza- 
tion forms and their spatial limits, which cannot be constrained by the existing data. Taking 
the CIE radiative cooling function ~ 5 x 10~ 23 ergs cm 3 s _1 at T ~ 2.4 x 10 6 K (Suther- 
land & Dopita 1993), we obtain the corresponding cooling rates as ~ 5 x 10 40 ergs s _1 and 
~ 7 x 10 39 ergs s -1 , which, within their uncertainties, are comparable to those obtained from 
observations of nearby disk galaxies of similar sizes (e.g., Wang et al. 2003). The cooling 
probably accounts for only a small fraction of the expected supernova (SN) mechanical en- 
ergy input in the Galaxy (~ 10 42 ergs s _1 , assuming 10 51 ergs per SN and one SN per 30 
years) . 

A more realistic spatial distribution of the hot gas may be between the above two ex- 
treme characterizations. Also our model fits are probably quite biased because of the limited 
number of the sight-lines. Nevertheless, the fitted values of z h and R h consistently indicate 
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D 



Rh 
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Fig. 9. — Left panel: The total hot gas column density as a function of the off-GC impact 
radius R p in the spherical distribution model. Njj is calculated via an integration along the 
dotted line (i.e., r = — oo to oo) in Fig. 8. The column density measurements and their errors 
have been enlarged by the ratio of to the model N H between the Sun and the source 
(see Fig. 8). This plot is for demonstration only; the actual fit is based on the comparison 
between the measured and modeled Nh values. The meaning of the symbols is the same as 
that in Fig. 7. The dashed line represents the best fit function with n c = 6.2 x 1CT 3 cm -3 , 
R h = 2.3 kpc (Eq. 16). Right panel: R h vs. n c contour at the 68%, 90%, and 99% confidence 
levels. 
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a small scale of the hot gas distribution. Because the fits incorporate the measurements to- 
ward the two extragalactic sources, this small scale suggests that a bulk of the highly-ionized 
X-ray-absorbing gas arises from regions close to the Galactic disk. 

In this work, we adopt the element abundances from Wilms et al. (2000) in which the 
abundance ratios of O/H and Ne/H are 0.58 and 0.71 times solar values given by Anders & 
Grevesse (1989). The real ISM abundances are still very uncertain (e.g., Savage & Sembach 
1996). However, if the real abundances are smaller or larger by a factor than those we 
adopted, the derived Njj (in §4), no, n c , and the total mass of the hot gas will be larger or 
smaller by the same factor accordingly, while the inferred T (in §4), Zh, and Rh are still valid. 
Since the Ne IX and the O VII/O VIII absorption lines are used here as the diagnostic of 
the ionization states and the gas temperature (Table 4), the results depend on the assumed 
Ne/O abundance ratio. However, this dependence can be partially compensated by a change 
in the best-fit T (please refer to Fig. 1). Therefore, the N H measurements should not be very 
sensitive to the assumed Ne/O ratio. In general, both T and Ne/O could be constrained 
simultaneously if multiple high S/N absorption lines of these two elements are available (§1). 



7. Comparisons with O VI Absorption and X-ray Emission Measurements 

We compare the above results from the X-ray absorption lines with those from obser- 
vations of Galactic far-UV O VI absorption lines and diffuse soft X-ray emission. Although 
the O VI absorption doublet (1031.93 and 1037.62 A) has commonly been detected in the 
spectra of both extragalactic and Galactic sources, the nature of the absorbing gas is still 
very uncertain. The O VI population peaks at T ~ 3 x 10 5 K, near the peak of the cooling 
curve of a CIE plasma. Therefore, such intermediate-temperature gas is expected primarily 
at the interfaces between the cool/warm clouds and the hot gas, consistent with the finding 
that about half of distinct high-velocity O VI absorbers have H I counterparts (Sembach 
et al. 2003). The bulk of the observed Galactic O VI density can be well described as a 
patchy exponential distribution with a scale height between 2.3 and 4 kpc (e.g., Zsargo et 
al. 2003). This scale height is greater than that inferred from our disk modeling of the X-ray 
absorption lines, which may indicate an evolutionary sequence of the hot gas. 

The X-ray-absorbing gas should also be responsible for some of the O VI population. 
We use the constrained n and the CIE ionization fraction of O VI at T ~ 2.4 x 10 6 K to 
estimate the O VI midplane density as 2.6(1.9, 3.5) xlO -9 cm -3 , which is ~ 1/7 of the value 
obtained from an O VI survey (Jenkins et al. 2001). Therefore, a significant fraction of 
the O VI absorption could arise from the X-ray-absorbing hot gas. Interestingly, the tightly 
constrained mean velocity dispersion of the O Vl-absorbing gas (Savage et al. 2003; Zsargo et 
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al. 2003) agrees well with the oxygen thermal velocity (~ 50 km s" 1 ) at T ~ 2.4 x 10 6 K, but 
is substantially greater than the velocity (~ 18 km s" 1 ) at the O VI peak temperature. Of 
course, much of the O VI line velocity dispersion may arise from turbulent and/or differential 
motion (e.g., due to the circular rotation of the disk and halo of the Galaxy). Detailed cross- 
correlation between the O VI and X-ray absorbing gas along individual sight-lines may 
provide insights into the connection between the O VI- and X-ray-absorbing gases. 

Fig. 2 shows a ROSAT all-sky diffuse 3/4- keV background map, including the locations 
of our studied sources and the two extragalactic sources. Among the seven LMXBs with the 
detected X-ray absorption lines (Table 4), five (X1820-303, X1728-169, X1636-536, X1735- 
444, and X1702-363) are within a 30° radius of the GC; in this Galactic inner region, the 
diffuse SXB is greatly enhanced in both 3/4 and 1.5 keV bands. The bulk of this enhancement 
may arise from the Galactic bulge (e.g., Wang 1998; Almy et al. 2000). The three sources 
(XI 118+480, X0614+019, and X0929-314), in which no absorption line is detected, are all 
toward outer Galactic regions with relatively low SXB emission; X1118+480 and X0614+019 
are also very close to us. Our estimated temperature and density of the X-ray-absorbing 
hot gas are between those inferred from the modeling of the X-ray emission gas in the 
solar neighborhood and the Galactic bulge (Wang 1998; Kuntz & Snowden 2000). These 
apparent correlation and consistency suggest that the absorption lines and the SXB emission 
arise from the same hot gas. 



8. Summary and Conclusions 

We have systematically searched for absorption lines produced by highly ionized species 
(O VII, O VIII, and Ne IX) in the spectra of 10 LMXBs observed with Chandra grating 
instruments. Much of our analysis is based on the multiplicative absorption line model, 
absline 2 , which we have constructed to allow for both an accurate absorption line profile 
fit and a joint analysis of multiple line transitions from same and/or different ions in a CIE 
plasma. 

We detect significant Ne IX K a absorption lines in seven of the 10 LMXBs. Three of 
these seven sources also show evidence for absorption by O VII K Q , Kg, and/or O VIII K a . 
The detected absorption lines most likely originate in the intervening diffuse hot gas, rather 
than in the LMXBs themselves, because of the absence of the expected line centroid shift 
and broadening as well as the lack of significant Ne X absorption and emission lines, which 
may be expected from a disk wind. We estimate the average temperature and the hydrogen 



2 The absline model, as implemented in XSPEC, can be obtained from the authors of this paper. 
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density of the absorbing gas to be ~ 2 x 10 6 K and ~ 4 x 10 3 cm 3 , under the assumptions 
of the CIE, unity filling factor, and ISM abundances. 

The hot gas is apparently located in and around the Galactic disk and preferentially 
in inner regions of the Galaxy consistent with the distribution of the diffuse SXB. We have 
included the column density measurements inferred from the absorption lines detected in 
two extragalactic sources to constrain the global distribution of the X-ray-absorbing hot 
gas. Modeled with an infinite disk with an exponential vertical distribution, the gas has a 
midplane density of ~ 6 x 1CT 3 cm -3 and a scale height of ~ 1 kpc. The gas, if indeed in a 
CIE state, can account for ~ 1/7 of the O VI absorption observed in the Galaxy. Modeled 
with a spheric radial distribution, the X-ray-absorbing hot gas has a predicted density of 
~ 6 x 1CT 3 cm" 3 at the GC and a core radius of ~ 2 kpc. But this spherical model gives a 
less satisfactory fit than the disk characterization. The small characteristic scales indicate 
that the bulk, if not all, of the cz ~ X-ray-absorbing hot gas detected in extragalactic 
sources to date is of Galactic origin. 

While far-UV absorption and X-ray emission studies have had a long history, X-ray 
absorption spectroscopy of the Galactic diffuse hot gas has just begun. Nevertheless, the 
results presented above demonstrate the potential of X-ray absorption line spectroscopy as 
a powerful diagnostic tool for probing the hot gas in our Galaxy. 

We thank D. Dewey for his detailed comments on an early version of the paper. We are 
also grateful to J. Kaastro, F. Paerels, T. Tripp, and T.-T. Fang for useful discussions and 
comments on the work, which is supported in part by NASA under grants AR4-5004A and 
G04-5046B. 
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